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Abstract 

Porcine deltacoronavirus (PDCoV) was first detected in Hong Kong and has recently spread to many countries around the 
world. PDCoV causes acute diarrhea and vomiting in pigs, resulting in significant economic losses in the global pork industry. 
In this study, a Chinese PDCoV strain, designated CHN-HG-2017, was isolated from feces of a suckling piglet with severe 
watery diarrhea on a farm located in central China. Subsequently, the virus was identified by an indirect immunofluorescence 
assay and electron microscopy. A nucleotide sequence alignment showed that the whole genome of CHN-HG-2017 is 97.6%- 
99.1% identical to other PDCoV strains. Analysis of potential recombination sites showed that CHN-HG-2017 is a possible 
recombinant originating from the strains CH/SXD1/2015 and Vietnam/HaNoi6/2015. Furthermore, the pathogenicity of this 
recombinant PDCoV strain was investigated in 5-day-old piglets by oral inoculation. The challenged piglets developed typical 
symptoms, such as vomiting, anorexia, diarrhea and lethargy, from 1 to 7 days post-inoculation (DPI). Viral shedding was 
detected in rectal swabs until 14 DPI in the challenged piglets. Interestingly, high titers of virus-neutralizing antibodies in 
sera were detected at 21 DPI. Tissues of small intestines from CHN-HG-2017-infected piglets at 4 DPI displayed significant 
macroscopic and microscopic lesions with clear viral antigen expression. Our analysis of the full genome sequence of a 
recombinant PDCoV and its virulence in suckling piglets might provide new insights into the pathogenesis of PDCoV and 
facilitate further investigation of this newly emerged pathogen. 


Introduction 


Porcine deltacoronavirus (PDCoV), a member of the fam- 
ily Coronaviridae, order Nidovirales, 1s an enveloped, sin- 
gle-stranded positive-sense RNA virus [1]. The genome of 
PDCoV 1s 25.4 kb in length, making it the smallest genome 
among the known coronaviruses (CoVs) [2]. The coding 
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region of PDCoV is flanked by short 5’ and 3’ untranslated 
regions (UTRs), and the gene order is 5’-UTR-ORF]la- 
ORF1b-S-E-M-NS6-N-NS7-3'-UTR. PDCoV encodes at 
least four structural proteins, including the spike (S), enve- 
lope (E), membrane (M), and nucleocapsid (N) proteins, as 
well as two accessory proteins, NS6 and NS7 [3-6]. The S 
protein is responsible for receptor binding and membrane 
fusion and acts as the major antigen inducing neutralizing 
antibodies [7]. The N protein is highly conserved and plays 
a critical role in viral RNA encapsidation [8]. The M and E 
proteins are important for virion assembly and budding [9]. 

During a molecular surveillance study performed by 
Yuen and coworkers in 2012, PDCoV was first identified 
in Swine specimens in Hong Kong [10]. Following the first 
detection of PDCoV in pigs with diarrhea in Ohio, USA, in 
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February 2014, dozens of PDCoV strains were identified in 
the United States, Canada, South Korea, China, and Thailand 
[4, 5, 11-16]. Piglets infected with PDCoV are characterized 
by watery diarrhea, variable vomiting, and dehydration [3]. 
Notably, PDCoV-positive piglets are frequently co-infected 
with other enteric viruses, such as porcine epidemic diar- 
rhea virus (PEDV), rotavirus, transmissible gastroenteritis 
coronavirus (TGEV), and porcine respiratory coronavirus 
(PRCV), which considerably increases the morbidity and 
mortality of piglets [17]. PDCoV outbreaks have been 
reported in many countries worldwide and have caused sig- 
nificant economic losses in the pork industry. Moreover, due 
to the unique genome replication strategy of coronavirus, 
PDCoV undergoes RNA mutation and recombination events 
at high frequency, giving rise to new variants, which makes 
it extremely difficult to control and eliminate [18]. 

In 2015, Saif and coworkers reported that a PDCoV strain 
named OH-FD22 caused severe atrophic enteritis accompa- 
nied by severe diarrhea and vomiting in gnotobiotic piglets 
[19]. Macroscopic examination revealed that the jeyunum 
and ileum are the major targets for PDCoV infection, and the 
pathological lesions were similar to, but milder than, those 
caused by PEDV [19]. Subsequently, the pathogenicity of 
several American strains, including USA/IL/2014, Michi- 
gan/8977, and Ohio CVM I, was investigated in piglets [20]. 
These American strains exhibited pathogenic characteris- 
tics similar to those of OH-FD22. In addition, it has been 
reported that newborn piglets inoculated with two recent 
Chinese PDCoV isolates, named CHN-HN-2014 and CHN- 
GD-2016, developed clear clinical signs, including severe 
diarrhea, vomiting, and dehydration [21, 22]. 

In this study, we successfully isolated a high-growth 
PDCoV strain, hereafter referred to as CHN-HG-2017, 
from feces of a pig on a farm located in Huanggang, Hubei 
province. Analysis of the genomic sequence suggested that 
CHN-HG-2017 is a potential recombinant virus derived 
from Chinese and Vietnam PDCoV strains. Furthermore, 
the pathogenicity of CHN-HG-2017 in 5-day-old pig- 
lets was investigated by clinical assessment, the quantifi- 
cation of viral shedding and distribution, histology, and 
immunohistochemistry. 


Materials and methods 
Collection of clinical samples and treatment 


Feces (n=17) from pigs experiencing diarrhea were 
obtained from farms in Huanggang, Hubei province. 
They tested PDCoV-positive by N-gene-based RT-PCR as 
reported previously [12]. These clinical samples were subse- 
quently subjected to virus isolation. Briefly, 1 g of PDCoV- 
positive feces was suspended in 5 ml of Dulbecco’s modified 
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Eagle’s medium (Gibco, Grand Island, NY, USA), vortexed 
for 5 minutes, and centrifuged at 4000xg for 15 minutes at 
4 °C. The supernatants were filtered through a 0.22-um filter 
(Millipore, Billerica, MA) and stored at -80 °C for PDCoV 
isolation. 


Virus isolation, propagation, and titration 


LLC-PK1 cells (ATCC CL-101) were obtained from the 
American Type Culture Collection (Rockville, MD, USA) 
and were used for PDCoV isolation. LLC-PK1 cells were 
cultured in DMEM supplemented with 10% fetal bovine 
serum (FBS, Gibco). The maintenance medium for PDCoV 
propagation was DMEM with 10% tryptose phosphate broth 
(TPB) and 10 pg of trypsin (Sigma) per ml. Cells were cul- 
tured in a 6-well plate until they reached 70-80% confluence 
and were then washed three times with PBS. Filtered inocu- 
lum (500 ul) in 1.5 ml of maintenance medium was added 
to the cell monolayers, and the cells were kept at 37 °C in 
5% CO, for 2 h. The cells were then washed three times 
with PBS, and 2 ml of maintenance medium was added to 
each well. The inoculated cells were cultured continuously at 
37 °Cin5% CO,. When a cytopathic effect (CPE) was obvi- 
ous in the infected cells, the plates were frozen and thawed 
twice. Subsequently, the supernatants were stored at -80 °C 
as seed stocks for plaque purification and for the next pas- 
sage. The isolates were plaque-purified as described previ- 
ously [20]. The virus titer was determined based on the 50% 
tissue culture infectious dose (TCID5,) on LLC-PK1 cells in 
96-well plates as described previously [20]. 


Immunofluorescence assay (IFA) 


The N gene of CHN-HG-2017 was amplified by RT-PCR, 
and the PCR product was cloned into the pET28a vector. 
The recombinant plasmid was expressed in Escherichia coli 
Rosetta (DE3) under IPTG induction. Subsequently, female 
BALB/c mice were immunized with the purified recombi- 
nant N protein. After hybridoma cell fusion and selection, 
one strain of hybridoma cells that secreted an anti-N protein 
monoclonal antibody (mAb), named A3, was selected. 
LLC-PK1 cells were seeded on sterilized coverslips 
placed in 24-well plates and then mock-infected or infected 
with plaque-purified PDCoV at a multiplicity of infection 
(MOI) of 0.01. At 24 h postinfection (HPI), cells were fixed 
with 4% paraformaldehyde for 15 min at room temperature 
(RT) and then permeabilized with 0.1% Triton X-100 for 
10 min. Subsequently, cells were blocked with 5% skimmed 
milk in PBS for 1 h at room temperature and incubated with 
mAb A3 for 1 h. The cells were treated with fluorescein iso- 
thiocyanate (FITC)-labeled goat anti-mouse IgG secondary 
antibodies (Invitrogen) and then stained with 4',6-diamidino- 
2-phenylindole (Invitrogen) for 15 min at room temperature. 
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Fluorescence was examined using a fluorescence microscope 
(Olympus [X73, Japan). 


Western blot analysis 


LLC-PKI1 cells were infected with PDCoV and harvested 
with lysis buffer (65 mM Tris-HCl] [pH 6.8], 4% sodium 
dodecyl sulfate, 3% DL-dithiothreitol and 40% glycerol) 
supplemented with PMSF. Proteins isolated from the lysate 
were separated by 12% SDS-PAGE and then transferred 
to a polyvinylidene difluoride membrane. The membrane 
was blocked with 10% dry milk and then incubated with the 
PDCoV-N-specific mAb A3. After washing three times with 
PBST, the membrane was incubated with horseradish per- 
oxidase (HRP)-conjugated goat anti-mouse IgG (ABclonal) 
for 45 min at room temperature. Signals were detected using 
a SuperSignal West Pico Luminal Kit (Pierce). 


Electron microscopy 


After inoculation with CHN-HG-2017, LLC-PK1 cells were 
harvested when more than 80% CPE was observed. Cells 
were frozen and thawed three times and then centrifuged at 
10,000xg at 4 °C for 20 min. The supernatants were filtered 
and ultracentrifuged at 30,000xg for 2 h at 4 °C. The pellets 
were resuspended in PBS and centrifuged in a non-linear 
30%-60% sucrose gradient at 30,000xg for 3 h at 4 °C. 
The virions that sedimented in the interphase between the 
30% and 45% sucrose layers were diluted with sterile PBS. 
The sucrose was then removed by centrifugation for 1 h at 
30,000Xg. The virions were resuspended in PBS, stained 
with 2% phosphotungstic acid, and examined using a trans- 
mission electron microscope (Hitachi H-7000FA, Japan). 


Phylogenetic and recombination analysis 


Viral RNA was extracted from CHN-HG-2017-infected 
LLC-PK1 cells using TRIzol Reagent (Invitrogen, USA) 
according to manufacturer’s instructions. cDNA was syn- 
thesized from the extracted RNA using random primers and 
a Takara RNA PCR Kit (Takara, Japan). The whole genome 
of CHN-HG-2017 was sequenced using 26 pairs of primers 
reported by Lorsirigool et al. [23]. Based on the sequence 
that was obtained, a pair of primers (PDCoV-S'GSP, 5'- GAT 
TACGCCAAGCTTTAGATCGCAAGGTGGGTT GGG 
CTCC-3'; PDCoV-3'GSP, 5’-GATTACGCCAAGCTTCGT 
CGTAAGACCCAACATCAAGCT-3’) was designed to 
verify the 5’- and 3’- terminal sequences of CHN-HG-2017 
using a SMARTer RACE 5’/3' Kit (Takara) according to 
manufacturer's instructions. The sequences were combined 
using DNASTAR Lasergene 7.0, and the assembled genome 
sequence was submitted to the GenBank database under the 
accession number MF095123. Multiple sequence alignments 


and phylogenetic analysis were performed using the neigh- 
bor-joiing method with MEGA 7 software (http://www. 
mehasoftware.net/). The aligned sequences were scanned 
for recombination events, using the Recombination Detec- 
tion Program (RDP). 


Inoculation of pigs with PDCoV strain CHN-HG-2017 


Animal experiments were carried out strictly according to 
the Guide for the Care and Use of Animals in Research of 
the People’s Republic of China. All experimental proce- 
dures were approved by the Scientific Ethics Committee of 
Huazhong Agricultural University (HZAUSW-2018-005). 
Twelve 5-day-old piglets were purchased from a commer- 
cial pig farm and were randomly divided into two groups, 
which were housed in separated containments. Prior to 
PDCoV inoculation, piglets were confirmed to be negative 
for PDCoV, PEDV, TGEV, and rotavirus by testing rectal 
swabs by RT-PCR. Neutralizing antibodies against PDCoV 
were not detected in sera of the piglets or sows. After accli- 
mation for 24 h, the challenged group was inoculated orally 
with a CHN-HG-2017 preparation with a titer of 1 x 10° 
TCID;,/ml (10 ml per pig), while the mock-infected group 
was inoculated with 10 ml of maintenance medium. Clinical 
symptoms, including vomiting, diarrhea, lethargy, and appe- 
tite loss, were observed and recorded daily. Rectal swabs 
were collected daily from each piglet from the first day post- 
inoculation (DPI) to the end of the experiment. Swabs were 
submerged into 1 ml of PBS immediately after collection. 
Randomly selected piglets from each group (n=3) were sub- 
jected to necropsy at 4 DPI, and prolonged PDCoV shedding 
was monitored in the remaining pigs until 21 DPI. At nec- 
ropsy, tissues of duodenum, jejunum, ileum, cecum, colon, 
rectum, liver, spleen, lung, kidney, and blood were collected 
for analysis by real-time RT-PCR. Tissues of duodenum, 
jejunum, ileum, cecum, and colon were also fixed with 10% 
formalin for histopathology and immunohistochemistry 
analysis. Sera collected at 0, 4, 7,21 DPI from piglets were 
tested for neutralizing antibodies against PDCoV. 


Real-time RT-PCR analysis 


The number of viral RNA copies in fecal swabs, sera, 
and tissue samples was determined by real-time RT-PCR. 
Briefly, fecal swab samples diluted in DMEM were homog- 
enized and centrifuged at 6000xg for 5 min. An aliquot 
(200 ul) of supernatant was extracted for viral RNA iso- 
lation. One gram of various tissue samples was mechani- 
cally homogenized in 1 ml of DMEM, and 200 ul of the 
supernatant was used for RNA extraction. Two hundred ul 
of serum was also used for RNA extraction. Viral RNA was 
extracted using TRIzol Reagent (Invitrogen, USA) according 
to the manufacturer's protocols and was treated with DNase 
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I. One pg of total RNA was used for cDNA synthesis by 
reverse transcription using an iScript cDNA Synthesis Kit 
(Bio-Rad). cDNA amplification was performed using the 
1Taq Universal Probes Supermix (Bio-Rad) with primers and 
probes targeting the N gene of PDCoV, as described by Ma 
et al. [3]. The recombinant plasmid pET28a-N containing 
the N gene of PDCoV was serially diluted tenfold to gener- 
ate a standard curve for each plate. Real-time RT-PCR was 
carried out using an Applied Biosystems 7500 Fast Real- 
Time PCR System (Life Technologies, USA). The thermal 
cycling parameters were as follows: 2 min at 50 °C, 10 min 
at 95 °C, and 40 cycles of 15 s at 95 °C and 1 min at 60 °C. 
The number of copies of PDCoV viral RNA in the tested 
samples was calculated based on the standard curve. 


Neutralization assay (NA) 


The presence of PDCoV-specific neutralizing antibodies in 
serum samples was determined using a virus neutralizing 
(VN) test. Briefly, LLC-PK1 cells were grown in 96-well 
plates until they reached 90% confluence. The CHN- 
HG-2017 stock was diluted in DMEM to make 200 TCID,;, 
viral stocks. Fifty pl of the diluted virus was then mixed with 
an equal volume of twofold serial dilutions of individual 
inactivated sera in 96-well plates and incubated at 37 °C 
for 1 h. After washing twice with maintenance medium, 
LLC-PKI1 cells were incubated with the mixture at 37 °C 
for 1 h. The mixture was then removed, and the cell mon- 
olayers were washed three times with PBS and cultured in 
maintenance medium at 37 °C in a 5% CO, incubator for 
2 days. The neutralization antibody titers were calculated 
as the reciprocal of the highest dilution of serum at which 
PDCoV infection was blocked. 


Histology and immunohistochemistry 


At necropsy, tissues of duodenum, jejunum, ileum, cecum, 
colon, and rectum collected from the inoculated and con- 
trol groups were separated and fixed in 10% formalin for 
24 h at room temperature and then dehydrated, embedded, 
sectioned, and mounted on glass slides. The samples were 
stained with hematoxylin and eosin (H&E), and the slides 
were examined and analyzed by conventional microscopy. 
Three representative villi and crypts with integrated longi- 
tudinal sections were randomly selected from each H&E- 
stained jejunum and ileum. Villus-height-to-crypt-depth 
(VH/CD) ratios were calculated using a computerized 
image system (NIS-Elements 3.2). Five-micron-thick sec- 
tions obtained from formalin-fixed paraffin-embedded tis- 
sues were placed onto positively charged glass slides, and 
the slides were air dried for 30 min. The tissue sections were 
deparaffinized, and the slides were rinsed and then incu- 
bated with target retrieval solution (Sigma-Aldrich, USA). 
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The sections were blocked with 1% BSA, incubated with 
PDCoV-N-specific monoclonal antibodies for 30 min, and 
then with HRP-conjugated goat anti-mouse IgG secondary 
antibodies (Zhongshan Biotechnology Co., Beying, China). 
The samples were visualized using a substrate-chromogen 
reagent (3-amino-9-ethylcarbazole/H,O,, Solarbio, China). 
Tissues from piglets in the negative control group were used 
as controls. 


Results 
Virus isolation 


Of the 17 PDCoV-positive feces inoculated onto LLC-PK1 
monolayers, only one caused obvious CPE in the form of 
enlarged, rounded and clustered cells at 24 HPI and rapidly 
detached from the monolayers at 36 HPI (Fig. 1A). This 
PDCoV isolate, designated as CHN-HG-2017, tested posi- 
tive for PDCoV and negative for other enteric viruses by 
RT-PCR (data not shown). CHN-HG-2017 was further pas- 
saged in LLC-PK1 cells for a total of 10 passages, indicat- 
ing that the PDCoV CHN-HG-2017 strain could be stably 
propagated in LLC-PK1 cells. 

To confirm that virus propagation had occurred, the 
infectivity of the plaque-purified CHN-HG-2017 strain in 
LLC-PK1 cells was tested by IFA and western blot with an 
mAb against the PDCoV N protein. As shown in Fig. 1A, 
PDCoV-specific immunofluorescence was detected in most 
cells at 24 HPI, and the N protein was mainly localized in the 
cytoplasm. N protein was also detected in PDCoV-infected 
LLC-PK1 cells by western blot (Fig. 1B). The morphology 
of the PDCoV particles purified from the supernatants of 
CHN-HG-2017-infected LLC-PK1 cells was examined by 
electron microscopy (EM) using negatively stained prepara- 
tions. As shown in Fig. 2, typical crown-shaped coronavi- 
rus particles with spiky surface projections and an average 
diameter of 80-160 nm were observed. These results clearly 
showed that a strain of PDCoV had been successfully iso- 
lated and identified. 


Phylogenetic and recombination analysis 
of the genome sequence of CHN-HG-2017 


The complete nucleotide sequence of the genome of CHN- 
HG-2017 was determined and deposited in the GenBank 
database under the accession number MFO095123. Similar 
to other PDCoV strains, the complete genome sequence 
of CHN-HG-2017 is 25,399 nt in length, consisting of 5’ 
UTR-1ab-S-E-M-NS6-N-NS7-3' UTR, and predicted tran- 
scription regulatory sequences (TRSs) were located at the 5’ 
end of each gene with a highly conserved core sequence of 
5’ -ACCAC-3’ (Table 1). Moreover, CHN-HG-2017 shared 
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Fig. 1 Cytopathic effect A 
(CPE) and IFA staining of 


LLC-PK1 cells infected with PDCovV- 
CHN-HG-2017. (A) Infected infected 
LLC-PK1 cells at 24 HPI (a) 

and 36 HPI (b) and mock- LLC-PKI1 
infected LLC-PK1 cells at 24 cells 

HPI (d) and 36 HPI (e) are 

shown. PDCoV-infected (c) 

or mock-infected LLC-PK1 

cells (f) were examined by IFA Mock- 

at 24 HPI using monoclonal . 

antibodies against PDCoV N infected 
protein. (B) PDCoV-infected LLC-PK1 
or mock-infected LLC-PK1 cells 

cells were analyzed at 24 HPI 

using a western blot assay with 

monoclonal antibodies against B 


the PDCoV N protein 


Fig. 2 Electron microscopic images of purified PDCoV particles. 
Bar, 50 nm 


97.6%-99.1% nucleotide sequence identity with the other 
74 PDCoV strains in the database. Phylogenetic analysis 
indicated that the PDCoV strains from the United States 
and South Korea clustered into a large clade, whereas CHN- 
HG-2017 clustered with other PDCoV strains detected in 
China since 2014. These data imply that the Chinese strains 
might share a more recent common ancestor with the US 
and South Korean strains than the ancestor of the strains 
from Thailand, Laos, and Vietnam (Fig. 3A). Notably, 


phylogenetic trees based on the S gene displayed a similar 
clustering pattern (Fig. 3B). 

The RNA genome of coronaviruses can undergo high- 
frequency recombination due to its large size [24]. We there- 
fore analyzed the genome sequence of CHN-HG-2017 using 
RDP. As depicted in Fig. 4A, two PDCoV strains, Vietnam/ 
HaNoi6/2015 (GenBank accession number: KX834351) and 
CH/SXD1/2015 (GenBank accession number: KT021234), 
might represent the parent lineages of CHN-HG-2017. Simi- 
larity plot and bootscanning analyses predicted the potential 
breakpoints in CHN-HG-2017 to be at nt 3153 and 4249. We 
therefore divided the aligned sequences into three separate 
segments: one from nt | up to the breakpoint at position 
3153, another from the breakpoint at position 3153 to the 
breakpoint at position 4249, and the last one from the break- 
point at position 4249 to the end. A phylogenetic analysis 
of these segments using the maximum-likelihood algorithm 
suggested the CH/SXD1/2015 strain from China to be its 
major parent and the Vietnam/HaNoi6/2015 strain from 
Vietnam to be its minor parent (Fig. 4B), providing strong 
evidence that recombination events had occurred. 


Clinical observations in pigs challenged 
with CHN-HG-2017 


To investigate the pathogenicity of CHN-HG-2017, 5-day- 
old piglets were challenged with CHN-HG-2017 via oral 
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Table 1 Coding potential and ORF Location (nt) Length (nt) 
putative transcription regulatory 
sequences of the PDCoV- 
CHN-HG-2017 genome 
lab = 5540-19333 18794 
S 19315-22794 3480 


E 22788-23039 202 
M 23032-23685 654 
NS6 = 23685-23969 285 
N 23990-25018 1029 
NS7 = 24084-24686 603 


feeding (Table 2). All piglets from the negative con- 
trol groups were in good health throughout the course of 
study with no observed clinical signs (Fig. 5A and B). As 
expected, three of the six piglets in the challenged group 
showed diarrhea at 1 DPI, and more piglets developed vary- 
ing degrees of diarrhea, vomiting, lethargy and anorexia at 
2-6 DPI (Fig. 5C and D). Lethargy, anorexia and the progres- 
sion to watery diarrhea were most severe at 3-6 DPI, and the 
piglets gradually recovered thereafter (Table 3). Notably, no 
PDCoV-associated deaths were observed in this study. These 
results demonstrated that CHN-HG-2017 is pathogenic to 
newborn piglets. 


Fecal shedding, virus distribution, and neutralizing 
antibody against CHN-HG-2017 


To determine the duration of fecal shedding and to moni- 
tor the changing load of CHN-HG-2017, we attempted to 
detect viral RNA in rectal swab samples from | to 21 DPI by 
quantitative real-time RT-PCR. The viral RNA was detected 
In S1x Out of six rectal swab samples collected from the chal- 
lenged group from 1 DPI to 14 DPI. The PDCoV RNA copy 
number in homogenized rectal swab supernatants reached a 
peak of over 10° copies/ml at 2-4 DPI (Fig. 6A). Notably, 
no PDCoV RNA was detected in samples from piglets in the 
negative control throughout the study. 

The distribution of CHN-HG-2017 in various tissues 
was also examined by quantitative real-time RT-PCR in 
three pigs from each group sacrificed at 4 DPI. Viral RNA 
was detected in two of three duodenums (average 10°*° 
copies/g), three of three jejunums (average 10°”° copies/g), 
three of three ileums (average 10°* copies/g), three of three 
ceca (average 10°** copies/g), three of three colons (aver- 
age 10’”° copies/g), and three of three rectums (average 
10°”? copies/g). It was also detected in one of three liv- 
ers (10°-”° copies/g), one of three spleens (10°! copies/g), 
one of three kidneys (10°” copies/g) and one of three lungs 
(10° copies/g), but no viral RNA was detected in the blood 
(Fig. 6B). 
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Length (aa) Frame ‘TRS location (nt) TRS sequence (s) 
(distance in bases to 
AUG) 
6265 +1,4+1 75 ACACCA(459)AUG 
1160 +1 19169 ACACCA(140)AUG 
84 +1 22765 ACACCA(18)AUG 
218 +1 23006 ACACCA(20)AUG 
95 +2 23633 ACACCA(46)AUG 
343 +2 23977 ACACCA(7)AUG 
201 +1 23996 ACACCA(82)AUG 


To investigate when PDCoV neutralizing antibodies 
were detectable in the infected piglets, serum samples 
collected at 0, 4, 7, 21 DPI were tested for PDCoV-neu- 
tralizing antibodies. Virus-neutralizing (VN) antibodies 
in sera were detected as early as 4 DPI in sera of PDCoV- 
inoculated piglets, and the titers increased gradually and 
peaked at 21 DPI (Table 4). PDCoV-specific neutraliz- 
ing antibodies were not observed in sera from the control 
piglets. 


Gross pathology, histopathology, 
and immunohistochemistry 


To determine the gross pathological and histological changes 
in piglets infected with CHN-HG-2017, three of the six 
pigs in each group were selected randomly for necropsy at 
4 DPI. The small intestines, where yellow watery contents 
had accumulated, were clearly transparent, thin-walled, and 
gas-distended in the PDCoV-challenged pigs (Fig. 7A). No 
other lesions were observed in any other organs of the chal- 
lenged groups or in any of the organs of the negative control 
pigs (Fig. 7B). These results showed that the small intestine 
is the major target of CHN-HG-2017 infection. Microscopic 
lesions in the small intestine were also observed by micros- 
copy. Compared to those of the negative control piglets, the 
intestinal villi of the PDCoV-infected piglets were blunted 
and fragmented and exhibited atrophy and vacuolation, espe- 
cially in the jejunum and ileum (Fig. 7C and D and Table 5). 
Moreover, aggregation of lymphocytes and hyperemia in 
the intestinal lamina propria were also observed (Fig. 7C 
and D), but the intestines of the negative control were nor- 
mal (Fig. 7E and F). Consistent with the histopathological 
results, PDCoV antigen was detected by immunohistochem1- 
cal staining in the cytoplasm of the infected villous entero- 
cytes (Fig. 7G and H). No PDCoV antigen was detected in 
the negative controls (Fig. 7I and J). These results show 
that CHN-HG-2017 infection can cause intestinal lesions in 
newborn piglets, especially in the jeyunum and ileum. 
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Discussion 


PDCoV is a newly emerged enteropathogenic coronavirus 
that causes diarrhea, dehydration and death in neonatal 
piglets. Since PDCoV was first reported in Hong Kong, 
it has spread to many countries around the world, lead- 
ing to significant economic losses in the pork industry. 
In mainland China, the prevalence of PDCoV in pigs was 
reported to be 33.71% in Jiangxi province [17], 23.4% in 
Guangdong province [15], and 10.2% in Hebei province 
[25]. Although the prevalence of PDCoV in China has 
been widely reported, only a few PDCoV strains could be 
successfully isolated in cultured cells until now. 

In the present study, feces from PDCoV-infected pigs on 
farms located in central China were used to isolate PDCoV 
in LLC-PK1 cells. Typical CPE was observed in infected 
LLC-PK1 cells at 24 h after inoculation. Following plaque 
purification, the PDCoV strain was confirmed by IFA and 
electron microscopy. Notably, only one PDCOV strain was 
isolated from a total of 17 PDCoV-positive feces samples, 
indicating that the success rate for virus isolation was still 
relatively low. This low rate of PDCoV isolation might be 
attributed to a small amount of infectious virus in the fecal 
samples collected for virus isolation. In addition, the time 
at which the sample is collected might also contribute to 
the success of PDCoV isolation. 

To characterize the virus isolate, the complete genome 
of CHN-HG-2017 was sequenced and analyzed. While the 
genome organization is similar to those of other reported 
PDCoVs, the genome size of 25,399 nt is somewhat dif- 
ferent from the others. Moreover, we identified predicted 
TRSs at the 5’ end of each gene of CHN-HG-2017, each 
with a highly conserved core sequence of 5’-ACCAC- 
3', which appears to be unique to members of the genus 
Deltacoronavirus. Multiple sequence alignment and phy- 
logenetic analysis showed that all of the known PDCoV 
strains share a high level of nucleotide sequence identities. 
However, CHN-HG-2017 was most closely related to other 
Chinese PDCoV strains. Our phylogenetic analysis sug- 
gests that PDCoV strains from China, the United States 
and South Korea are clustered in the same clade, whereas 
the PDCoV strains from Thailand, Laos, and Vietnam are 
clustered in another clade. 

Like other coronaviruses, PDCoV has the potential to 
undergo recombination at a very high frequency. In this 
study, we employed the Recombination Detection Program 
(RDP) to analyze the genome of CHN-HG-2017 in search 
of any recombined fragments. Intriguingly, CHN-HG-2017 
was found to be a possible recombinant virus derived from 
parental strains from China and Vietnam. The major par- 
ent and minor parent were identified as strain SXD1 from 
China and strain HaNoi6 from Vietnam, respectively. The 


potential breakpoints in CHN-HG-2017 were located at 
nt 3153 and 4249 in the nsp3 region. Deletion mutations 
in the nsp3 region in PDCoV have been reported recently 
[6, 26], suggesting that nsp3 which encodes a papain-like 
protease, may be a hotspot for genetic modification. To our 
knowledge, this is the first time that a recombination event 
has been identified in the Nsp3 region of PDCoV. The 
emergence of recombinant strains from China and Vietnam 
will probably bring new challenges to the prevention and 
control of PDCoV. 

To investigate the pathogenicity of this putative recom- 
binant strain, 5-day-old piglets were inoculated orally with 
CHN-HG-2017. As expected, severe diarrhea with vomit- 
ing and lethargy in piglets was observed from 1 to 6 DPI, 
suggesting that CHN-HG-2017 is pathogenic to newborn 
piglets. The viral RNA distribution in various tissues was 
also examined. Viral RNA levels were higher in intestines 
than in other tissues. Previously, it was reported that coro- 
navirus HKUI15 may be able to spread through the respira- 
tory route, in addition to fecal-oral transmission [27]. In line 
with other reports, viral RNA of CHN-HG-2017 was also 
detected in the lungs of some animals, albeit with no gross 
lesions observed, suggesting that PDCoV may also cause 
respiratory infection. 

It is noteworthy that we only observed microscopic 
lesions in the jejunum and ileum in piglets challenged with 
CHN-HG-2017. In striking contrast, microscopic lesions 
were observed in all sections of small intestine infected with 
PEDV [28]. Moreover, unlike PEDV, PDCoV infection usu- 
ally does not result in the death of newborn piglets despite 
the fact that these animals exhibit lethargy and anorexia and 
have watery diarrhea. These results indicated that the viru- 
lence of PDCoV is milder than that of PEDV. 

Recently, it was shown that VN antibodies could be 
detected in 14-day-old piglets inoculated with PDCoV 
OH-FD22 at about 7 DPI and that VN antibody titers 
increased gradually and remained high at the end of the 
experiment [29]. However, to our knowledge, there are cur- 
rently no reports on the development of PDCoV-specific 
antibodies in sera of pigs infected with Chinese PDCoV. 
Our study showed that serum PDCoV VN antibodies could 
be detected as early as 4 DPI in PDCoV-infected piglets. 
VN antibody titers increased and remained high until 21 
DPI, when all pigs in the challenged groups had completely 
recuperated. 

In conclusion, we successfully isolated, characterized and 
obtained a high titer of CHN-HG-2017 in cultured cells. We 
showed that CHN-HG-2017 is probably a recombinant virus 
derived from CH/SXD1/2015 and Vietnam/HaNoi6/2015. 
Pathogenicity testing of CHN-HG-2017 revealed that this 
virus is still pathogenic in piglets despite multiple passages 
in cultured cells. High levels of serum VN antibodies were 
also detected in sera of infected animals. Not only will our 
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<Fig.3 Phylogenetic analysis using nucleotide sequences of the com- 
plete genome (A) and the S gene (B) of PDCoV strains from the 
GenBank database. CHN-HG-2017 is indicated by a triangle. The 
phylogenetic tree was constructed by the neighbor-joining method in 
MEGA 7 (http://www.megasoftware.net). A bootstrap analysis was 
performed with 1000 replicates, and the bootstrap values are indi- 
cated on each branch. The scale bar indicates nucleotide substitutions 
per site. The branches in pink, blue, yellow and green represent a 
clade from the United States and South Korea; a clade from China; 
a clade from Thailand, Laos, and Vietnam; and the outgroup, respec- 
tively 
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findings be useful for understanding the molecular epide- 
miology and pathogenesis of PDCoV, but they also provide 
insights for the development of effective vaccines against 
this important pathogen. 


95.0 
94.5 3153 4249 KX834351.1 
94.0 
0 2,000 4,000 6,000 8,000 10,000 12,000 14,000 16,000 18,000 20,000 22,000 24,000 26,000 ~_—-28,000 
Genomic positions nt 
- KP757891.1 
KX834351.1 
* MF095123.1 
Recombinant 3153 4249 
points nt 
Phylogenetic analysis: 1-3153 nt 3153-4249 nt 4250-end 
53) KJ462462.1 KU981061.1 
75; 3Q065043.2 MF431743.1 


JQ065048.1 


4 Potential recombinant MF095123.1; © Potential major parent KP757891.1; m Potential minor parent KX834351.1. 


Fig.4 Recombination analysis of the complete genome sequence 
of PDCoV strain CHN-HG-2017. (A) Possible recombination 
events were identified using the Recombination Detection Program 
(RDP). Potential recombination events supported by strong P-values 
(<107!°) were examined further using similarity plot and bootscan 


analyses implemented in Simplot 3.5.1. (B) Phylogenetic trees based 
on three different segments of the genome were constructed using the 
maximum-likelihood algorithm with the LG model with bootstrap 
values determined by 1000 replicates in the PhyML (version 3.0) 
software package 
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Table 2 Study design for oral 


. Group No. of piglets Inoculum No. of piglets used 
inoculation ok five-day-old inoculated for necropsy at 4 
piglets with PDCoV isolate DPI 
CHN-HG-2017 

PDCoV 6 PDCoV 10°TCID50/ml; 10 ml 3 

Negative control 6 Maintenance medium; 10 ml 3 


Fig.5 Clinical assessment of piglets challenged with CHN-HG-2017. Five-day-old pigs were examined 4 days after inoculation with PDCoV 
strain CHN-HG-2017 (C, D) or DMEM medium (A, B) 
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Table 3 Clinical observations in piglets challenged with PDCoV 


Watery 
diarrhea 


0/6 
1/6 
2/6 
4/6 
2/3 
1/3 
1/3 
0/3 
0/3 


CHN-HG-2017 
DPI Lethargy & Fecal consistency 
anorexia 
Normal Mild diarrhea 
0 0/6 6/6 0/6 
1 2/6 3/6 2/6 
Bi 4/6 2/6 2/6 
3 4/6 1/6 1/6 
4 3/3 0/3 1/3 
5 3/3 0/3 2/3 
6 2/3 1/3 1/3 
y 2/3 1/3 2/3 
8-21 0/3 3/3 0/3 
A Fecal shedding of 5-day-old pigs 
inoculated with PDCoV 
10 
= 9 
= 
g 8 
> § 
ax 5 
a 4 
oO 3 
D 2 
o 
— 1 
0 


Day post-inoculation 


Table 4 Anti-PDCoV neutralizing antibody titers in the serum of 
5-day-old pigs inoculated with PDCoV CHN-HG-2017 


Group Pig no. 
PDCoV C2 
C4 
G5 
Negative control NCI 
NC4 
NCS5 
B 


O DPI 


<1:4 
<1:4 
<1:4 
<1:4 
<1:4 
<1:4 


4 DPI 


7 DPI 21 DPI 


64 
64 
64 
<1:4 
<1:4 
<1:4 


inoculated with PDCoV at DPI 4 


3/3 3/3 «3/3 


Log10 RNA copies/g 
Oem NnNWAADN WO OC 


3/3 3/3 


256 
128 
128 
<1:4 
<1:4 
<1:4 


Viral RNA distribution of 5-day-old pigs 


Fig.6 Fecal viral shedding and virus distribution in PDCoV-challenged pigs. (A) Fecal viral shedding in pigs challenged with CHN-HG-2017. 


(B) Virus distribution at 4 DPI in pigs challenged with CHN-HG-2017 
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Fig. 7 Intestinal lesions in a PDCoV-challenged pig. (A) Macro- tissue section of a negative control pig. (G, H) Immunohistochemi- 
scopic lesions of a PDCoV-challenged piglet at 4 DPI. (B) Mac- cally stained jejunum (G) and ileum(H) tissue section of a PDCoV 
roscopic picture of a negative control piglet at 4 DPI. (C, D) challenged piglet. (I, J) Immunohistochemically stained jejunum (I) 
H&E-stained jejunum (C) and ileum (D) tissue section of a PDCoV- and ileum (J) tissue section of a negative control pig. Scale bars are 
challenged piglet. (E, F) H&E-stained jejunum (E) and ileum (F) shown in each picture 
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coronavirus-like particles by co-expression of viral envelope pro- 
tein genes. EMBO J 15:2020—2028 
10. Woo PC, Lau SK, Lam CS, Lau CC, Tsang AK, et al. (2012) 


Table 5 Histology and immunohistochemistry analysis of intestines 
infected with PDCoV strain CHN-HG-2017 


Group Pig no. WVH:CD, Antigen detection in the Discovery of seven novel mammalian and avian coronaviruses in 
Mean intestine Deltacoronavirus supports bat coronaviruses as the gene source of 
(+SD) —_-___—— Alphacoronavirus and Betacoronavirus and avian coronaviruses 
Duodenum Jejunum [leum as the gene source of Gammacoronavirus and Deltacoronavirus. 
JVI 06540-06511 

PDCoV 2 2.04 (0.5) + naa pees 11. Marthaler D, Raymond L, Jiang Y, Collins J, Rossow K et al (2014) 
3 2.3 (0.4) + ++ +++ Rapid detection, complete genome sequencing, and phylogenetic 

Negative control 2 5.5 (0.4) ; . analysis of porcine deltacoronavirus. Emerg Infect Dis 20:1347 
5 5.32 (0.9) - : ; 12. Wang L, Byrum B, Zhang Y (2014) Detection and genetic char- 


acterization of deltacoronavirus in pigs, Ohio, USA, 2014. Emerg 
Infect Dis 20:1227 

13. Wang L, Byrum B, Zhang Y (2014) Porcine coronavirus HKU15 
detected in 9 US states, 2014. Emerg Infect Dis 20:1594—1595 

14. Wang L, Zhang Y, Byrum B (2014) Complete genome sequence 
of porcine coronavirus HKU15 strain IN2847 from the United 
States. Genome Announc 2:e00291-00214 

15. Chen F, Zhu Y, Wu M, Ku X, Yao L et al (2015) Full-length 
genome characterization of Chinese porcine deltacoronavirus 
strain CH/SXD1/2015. Genome Announc 3:e01284-01215 

16. Janetanakit T, Lumyai M, Bunpapong N, Boonyapisitsopa S, Chai- 
yawong S et al (2016) Porcine deltacoronavirus, Thailand, 2015. 
Emerg Infect Dis 22:757 

17. Song D, Zhou X, Peng Q, Chen Y, Zhang F et al (2015) Newly 

emerged porcine deltacoronavirus associated with diarrhoea in 

Swine in China: identification, prevalence and full-length genome 

sequence analysis. Transbound Emerg Dis 62:575—580 
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